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We present a new approach to static and dynamical properties of holes and spins in weakly doped 
antiferromagnets in two dimensions. The calculations are based on a recently introduced cumulant 
approach to ground-state properties of correlated electronic systems. The present method allows 
to evaluate hole and spin-wave dispersion relations by considering hole or spin excitations of the 
ground state. Usually, these dispersions are found from time-dependent correlation functions. To 
demonstrate the ability of the approach we first derive the dispersion relation for the lowest single 
hole excitation at half-filling. However, the main purpose of this paper is to focus on the mutual 
• influence of mobile holes and spin waves in the weakly doped system. It is shown that low-energy 

spin excitations strongly admix to the ground-state. The coupling of spin waves and holes leads 
to a strong suppression of the staggered magnetization which can not be explained by a simple 
rigid-band picture for the hole quasiparticles. Also the experimentally observed doping dependence 
O I of the spin-wave excitation energies can be understood within our formalism. 
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^ . I. INTRODUCTION 

oo i 

The physical properties of high-temperature superconductors are strongly influenced by electronic correlations. The 
investigation of strongly correlated electronic systems has therefore become a major topic for the understanding of 
the high-T c .materials. The undoped compounds arc antifcrromagnetic Mott-Hubbard insulators. Neutron scattering 
experimentsM show the existence of spin-wave excitations which can be described by conventional spin-wave theory 
for the isotropic spin S = \ Heisenberg model on a square latticeB. The doped materials show a strong dependence 
of the magnetic properties on the hole concentration S in the CuC>2 planes. With increasing hole concentration both 
the Neel temperature and the staggered magnetization decrease and vanish at a critical hole concentration S c of a few 
percent before the system becomes paramagnetic and metallic (or superconducting at sufficiently low temperatures). 
Upon doping also the spin- wave velocity decreases and vanishes at approximately the same hole concentration 8 C . At 
' ^ , the same time the long-wavelength spin-wave modes become overdamped. 

The essential aspects of the low-energetic electronic, degrees of freedom of the CuC>2 planes are by now believed to 
q \ be well described by the two-dimensional t-J modelerE!: 

H = -* E (C^ + sW) + J E ^ - ! r 1 ) • « 

(ij)o- (ij) 



' Here, Si is the electronic spin operator and rii the electron number operator at site i. The symbol (ij) refers to a 
summation over pairs of nearest neighbors. In the following we denote the two antiferromagnetic sublattices by f and 
L In all sums we shall use i G| and j Gj.. Note that the Hamiltonian ([[]) is defined in the subspace of the unitary 
space without double occupations of sites. The electronic creation operators c ia are not usual fermion operators but 
rather exclude double occupancies: 

4 = 4(1-^). (2) 

At half-filling the t-J Hamiltonian reduces to the antiferromagnetic Heisenberg model. 



Experiments show that already a few percent of additional holes (away from half-filling) destroy the long-range 
antiferromagnetic order in the CuC>2 planes. This demonstrates the importance of the interplay between antiferro- 
magnetism and the motion of holes in the high-T c superconductors. To investigate the mutual influence one may 
start from the t-J Hamiltonian (|l|). However, due to strong correlations usual diagrammatic techniques based on 
Wick's theorem can not be easily applied to (|j). Neither the first nor the second part of (|l|) is bilinear in fermion 
operators, and the creation and annihiliation operators ct and Ci„ do not obey simple anticommutation relations. For 
this reason, non-standard analytical methods like variational wavefunctions, coupled-cluster methods, or slave-boson 
and slave-fermion techniqucsEnl3 have been employed. 
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In the following, we present a static approach to evaluate static and dynamical properties in weakly doped antifer- 
romagnets. The calculations are based on a cumulant method for computing the ground-state energy of correlated 
electronic systems. This is in contrast to the usual approach to dynamical properties which is based on dynamical 
quantities like time- or frequency-dependent correlation functions. Thejpaper is organized as follows: In Sec. II we 
shall describe the cumulant method which was recently proposed in refsBcl. In Sec. Ill this formalism is applied to the 
motion of a single hole generated in the groyuod .state of a quantum antiferromagnet at half-filling. This problem was 
already investigated by a number of authorsoll3't£rEil. Our ficst-aim is to show that our approach is able to reproduce 
results known from literature, see for instance review articlestj'lij. Processes leading to the hole motion are described 
within the concept of path operators which leads to the well-known spin-bag picture. For the case of one hole we 
obtain a quasiparticle dispersion which has minima at (±7r/2, ±7r/2) and a bandwidth of 1.4J...1.5J in agreement 
with several analytical and numerical calculations. Our main aim is to investigate the coupling of spins and holes in 
the weakly doped regime. This is the subject of Sections IV and V. In Sec. IV we describe the ground state for the 
case of small doping. Then we evaluate the staggered magnetization from the ground-state energy by introducing an 
external field coupling to the staggered magnetization. In Sec. V we derive the dispersion relation for spin waves in 
the doped system. We find that the spin wave energies are strongly renormalized due to the presence of holes. The 
spin- wave velocity vanishes at a critical hole concentration of a few percent. The strong coupling between spin waves 
and holes also explains the experimentally observed fast decrease of the staggered magnetization with increasing hole 
doping. 



II. CUMULANT APPROACH 



Conventional treatments of systems with electronic correlations start from the uncorrelated limit and include many- 
body effects by perturbation theory. This includes summation over classes of Feynman diagrams. For a treatment 
of strongly correlated electrons one would like to proceed the other way round, i.e., by starting from a local picture 
including the electronic correlations and by expanding with respect to the hybridization interactions. In this case 
the unperturbed Hamiltonian Hq contains two-particle operators. This leads to the difficulty that Wick's theorem 
is no longer applicable. Therefore, usual diagrammatic techniques cannot be used for strongly correlated electronic 
systems. For this reasomiither methods for treating correlated electrons have been developed, e.g., slave-boson and 
slave-fermion techniquesErEx n n 

An alternative approach for calculating expectation values and dynamical correlation functionsclcl is based on 
the introduction of cumulants. Provided that the Hamiltonian of the system can be split into H = Hq + Hi with 
eigenstates and eigenvalues of Hq known, this method uses the decomposition 

e -A* = e -X(H 1+Lo ) £ -XH (3) 

which can be proven by comparing the equations of motion of either side with respect to A. The Liouville operator 
Lq is a superoperator defined by LqA = [Hq,A]-. Let us denote the ground state of the unperturbed Hamiltonian 
Hq by |0 O ) and its energy by e 

Hq\4>q) =eo\<h). (4) 

For the following it is useful to introduce a cumulant bilinear form defined by 

(A\B) = {<f> \rfB\fo) : = {falAiBlfa)*, (5) 

where (...) c denotes a cumulant expectation value. For a detailed discussion of cumulants see e.g. KubcB 
Our aim is to calculate the ground state energy E of H: 

H\^q) = Eo\ih). (6) 

Using (||) one can show§ that Eq is given by 

E =e + (H 1 \n) = (H\Q.), (7) 

il = l+ lim -— -Hi. 

x^O x - {Lq + Hi) 

The operator f2 has similarity to the so-called wave operator (or Moeller operator known from scattering theory). It 
transforms the ground state |(/>o) of the unperturbed system into the exact ground state \ipo) of H. 
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The relation (0) can be applied to either weakly or strongly correlated systems because its use is independent ol 
the operator statistics , i.e., it is valid-,for fermions, bosons or spins. Based on the same approach alsg.-dynamical 
correlation functions can be calculated^. For a derivation and detailed discussion of relation see refs.Qfcl 

Treating cumulant expectation values one must distinguish between prime and composite operators. A prime 
operator is a single entity in the cumulant evaluation procedure. Expanding O given in (0) the resulting products of Lq 
and Hi are composite operators in the cumulant ordering. Assumingpthat f2 can be represented by an universal analytic 
function / of a prime operator S, i.e., £1 = f(S), it has been shownl!3 that must be of the form O = f(S) = expS*. 
As a generalization we shall use in the following an exponential form as an ansatz for f2, i.e., 

!] = e^ W " (8) 

with a finite set of relevant operators S u replacing S. The S v have to be chosen in such a way that exp(^ i/ X^SpM(f>o) 
(with appropriate parameters Aj,) represents a good approximation of the exact ground state. Following ref.113 we 
obtain a set of coupled equations for calculating the ground-state energy E : 

E Q = (H\e s ), 

0=(S fl \He s ), (9) 
S = S \ V S V . 

V 

The expansion coefficients X, can be determined from the equation ^ji. Note that the relation = (A\Hfl) holds 
for all operators A (see ref.E3). 

In several cases the set of nondinear equations (|^) can be transformed into an eigenvalue problem. First, the 
cumulant expectation values with the exponential can be evaluated according to appendix A. This leads to 

E = (He s ), 

(S+Hc s )=E (S+e s ), M =l,... (10) 

with S defined in (^) 3 . With the formal definition Sq = 1 and the assumption (£>„) = for all operators S v {y — 1, ...) 
in S we can formally include (|To|) i : 

(S+He s )=E (S+e s ), M = 0,l,... (11) 

If we now consider the case that any products of operators can be expressed by linear combinations of other 
operators from the set {S v } we can decompose the exponential into 

exp(£A„S„) =Y,P»S V . (12) 



where we have introduced a new set of parameters {13^}. Inserting (12) into (|ll|) leads to a generalized eigenvalue 
problem: 

Y,^{S+HS V ) = E 52Pu(S+S v ), m = 0,1,... . (13) 

t/=0 u=0 

Note that (|ll| ) is a priori size-consistent even if the sum in S is restricted to a finite set of operators S v due to the 
exponential which contains perturbations of the ground state up to infinite order. Note that, however, dT3| ) with a 
finite set {S„} is no longer size-consistent since the state (3 V S V \4>$) contains only a finite number of excitations. 
This is the same case as in a configuration-interaction calculation. 



III. HOLE MOTION (REVISITED) 

In this section we use the present formalism to describe one hole moving in an uud apftd-antifcrromagnet (for hole 
concentration S — > 0). This problem has already been studied in a number of papersBEjIljtj. The main contribution 
to the hole motion is the following: The motion of the hole locally destroys the antiferromagnetic spin order leading 
to a string of spin defects. However, quantum spin fluctuations can repair pairs of frustrated spins. This process 
leads to a coherent motion of the hole for each of the twp-|p*lplattices. Some of the above calculations are based on 
the motion in an ideal antiferromagnetic background, e.g.r 8 H 19 H 2 H. However, spin fluctuations in the antiferromagnetic 
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ground state allow for additional hole motion processes and thus influence the properties of the hole quasiparticlesc 
Several authors have also studied the hole motion using numerical methods like exact diagonalization and Quantum 
Monte Carlo methods. 

In the following wc discuss the hole dispersion relation by calculating the ground-state energy of a system at half- 
filling which contains one additional hole with fixed momentum k. Contrary to usual methods this approach is based 
on a static view-. of-,the hole motion problem. For a proper description of the one-hole states we use the concept of 
path operatoral3~LJ and define path concatenation operators A n = A n -\ + A n [ . The operators A n -\ and A n i refer to 
the two sublattices. A n operating on the Neel state with one hole, Ci^ij) n eel) , moves the hole n steps away and creates 
a path or string of n spin defects attached to the transferred hole. Explicitely, the operators A n ^ are defined by 



- 1 V- + r 
Al T = "7= '•■••'•<. /,> ' 



. V* ' 

M = . = = ^enSfc+RffRji, (14) 

V z o i z o ~ 1) ijlm 
( i et, je|, ieT, mei) 

The operators A n ^ for the 'down' sublattice are defined analogously with all spins reversed. z$= 4 denotes the number 

of nearest neighbor sites in the lattice. The matrices Rji and iJy allow the hole to jump to its four nearest neighbors 
in the first step and to only three new nearest neighbors by hopping forward in each further step: 

jjj _ f 1 i, j nearest neighbors , 
]l ~ 1 otherwise ^ ' 

_ f 1 j, I nearest neighbors and I 7^ % 
l i \ otherwise 

Let us split the Hamiltonian (jl]) into an unperturbed part H$ and into a perturbation H\ according to 
H = H Ismg = J (SiSj - T ^ L ) + J(N-2) , 

<ij> 

H x = H t + H ± (16) 

= — t 2^ i^Cjo- + Cj^Cia) + — ^ (S i Sj~ + S^Sj ). 

<ij>,cr <ij> 

The unperturbed Hamiltonian Hq is essentially the Ising part of the Heisenberg exchange in (|l|) whereas Hi contains 
the transverse part as well as the conditional hopping contribution. The ground state \4>o) of i^o with one hole with 
momentum k is given by 

I0o> = -7^7= £ eikH4c *Tl^rf> • ( 17 ) 

For the wave operator we use the exponential ansatz of Sec. II. In order to include the path operators A n described 
above and also to take into account ground-state spin fluctuations generated by H± we choose the following form 

00 

fl = cxp KA n ) exp (fiAp) (18) 

n=l 

Here, Ap is a spin-flip operator defined by 

A F = Y, s TSjPij (i €t, j el) , (19) 

<ij> 

P l j = JJ(nj-r +n n ). 
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It creates pairs of spin flips on nearest neighbor sites and corresponds to the well-known Bartkowski wavefunctionES 
which describes spin fluctuations in the ground state of the Heisenberg antiferromagnet alone. The projector 
prevents spin fluctuations next to the hole ( denotes all sites next to the pair of nearest neighbor spins ij ). It 
is introduced to avoid ambiguities since such an excitation is equivalent to a path of length 2 and is described by 
the path operator A^. Note that Ap does not commute with the path operators A n because of the presence of the 
projector Py. Note also that we have introduced in £1 (|lq) a product of two exponentials rather than one. This was 
done for simplification of further calculations and can be considered as an extension of (|^) . For the same reasons we 
shall also take into account spin fluctuations only up to the first order in Ap, i.e., 

exp(fxAp) w 1 + fj.Ap. (20) 

Using ( p0| ) we shall find — /i = | <C 1 which means that higher order terms in Ap are indeed neglegible. 

Following the method described in Sec. II we obtain the following set of non-linear equations to determine the 
ground-state energy Eq and the coefficients A„. 

oo 

E = (H\ cxp(J] X n A n ) (I + nA F )) , 

n=l 

oo 

0= ((l + fiA F )A l \H exp(J2*nA n )(l + fiA F ))), i = 1, 2, ... , (21) 

n=l 

oo 

= {Ap \H exp(J] X n A n ) (1 + fiAp)) . 
n=l 

Note that for algebraic reasons we have used composite operators A in the second of these equations (compare (0)). 
In (|2l]) all cumulant expectation values have to be taken with respect to the unperturbed ground state ([n]). As 
described at the end of Sec. II we use the relations given in appendix A to transform the set of coupled non-linear 
equations (21 ) 2 into a generalized eigenvalue problem. 

After evaluating the cumulants and expanding the exponential one obtains terms with products of path operators. 
For such a product of two path operators either the two operators couple to different holes creating two different 
paths, or the second path operator couples to the first one concatenating both paths. In the latter case both paths 
couple to the same hole. Therefore, their effect can be written as one path operator with both lengths added. Thus, 
after having evaluated the cumulants we introduce new path coefficients j3 n instead of the A„. Application of fl on 
the unperturbed ground state leads to 

oo oo 

exp(^A„A„) exp(nA F )\cl)o) = (1 + ^ /3 n A n )(l + ptA F )\<j) a ) (22) 

n—1 n—1 

There is a non-linear one-to-one correspondence between the (3 and A. In the following we only calculate the (3 n since 
the knowledge of the original A„ is not needed. 
From (^f|)i we find for the energy 

oo 

E = (H\J2 PnA n ) + f i(H\A F ) . (23) 

71=1 

Note that no mixed contributions in (3 n and /j, are obtained because H can only repair two spin defects created by Q. 
The first term in (|2^) depends on the hole momentum k and therefore describes the hole dynamics. The second term 
represents the ground state energy of a Heisenberg antiferromagnet (with one hole). The coefficients /i and (3i have 
to be determined from (pT|) 2.3- By evaluating the cumulants according to Appendix A, expanding the exponentials as 
described above and reintroducing new cumulants we find 

00 

( (I + tiA F )Ai\H(J2 PnA n ) (1 + fiAp) ) 

n=l 

00 00 

= {H\J20nA n ) x Hl+fiA F )Ai\(J2l3 n A n )(l+nA F )) (24) 

n—1 n—1 

OG OO 

+ ^(Ai\J2^A n ) x (A F \H{J2^A n )(l + ^A F )) 

n—1 n—1 
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and 



= (A p \eQ2 (3 n A n ) (1 + (iA F ) ) + (A F \H) + fi{A F \HA F ) . (25) 

n=l 

The set of equations ( p4| ) can be considered as a generalized eigenvalue problem for the (3 n and has to be-solved 
numerically (note (A n |j4 m ) = 6 nm ). The expectation values ace calculated analogously to a recent papered. All 
important contributions to the hole motion including spiral pathstJ and processes due to ground state spin fluctuations 
are taken into account. From ( |25|) we obtain /i = — i. ThifL value can also be derived directly from the expression 
for the wave operator f2 given in (Q) with Hi — H±, see ref£3. After having calculated the values for the coefficients 
j3i and \i we are able to evaluate the ground-state energy Eq of the system given by (|2^). The energy dispersion for 
different values of tj J is shown in Fig. 1. This dispersion relation agrees very well with other theoretical results and 
also with data found from exact diagonalization of small clusters. The energy minima are located at (±7r/2, ±7r/2), 
the total bandwidth is 1.4J...1.5J (see Fig. 2). ___ 

The path coefficients strongly decrease in space within a few lattice constants, see e.g.t3EifE3. The cloud of spin 
defects surrounding the hole (spin bag) is small with an average radius of approximately 2 lattice constants for tj J = 5. 
Therefore it is a reasonable approximation to include only the first path coefficients to describe the quasiparticle. 

The present approach can also be extended to finite hole concentrations. If we neglect hole-hole interaction processes 
we arrive at a rigid-band approximation. The Fermi surface shows the well-known hole pockets (see Fig. 3). 

Note that within the method presented here all quantities are calculated consistently. Especially, our results are 
not based on the calc ulati ons of the path coefficients for the localized hole casell3. This is an improvement compared 
to former calculationscSE3. There the hole dispersion relation problem was tackled in first order perturbation with 
respect to spin fluctuations with a zero order wave function, i.e., with the path coefficients of the localized hole. 



IV. GROUND STATE AND STAGGERED MAGNETIZATION OF THE DOPED SYSTEM 

Now we turn to the description of the spin dynamics in the doped system. We want to focus on the coupling 
between spin waves and hole motion. In the following S denotes the hole concentration away from half-filling. Our 
system with TV lattice sites possesses M = 8N dopant holes. 

The change of the magnetic jarouerties due to the presence of mobile holes was already investigated, within the 
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t — J model by several authorsEjEa. In some of the calculations the so-called magnetic polaron Jtiode^ 
used which can be derived from the t-J Hamiltonian ([!]) by; use of slave-fermion methods. In ref.E3 dynamical spin 
susceptibilities were calculated using projection techniquet^TEil instead. There a ground state according to a rigid- 
band approximation for the hole quasiparticles was used to evaluateLexpectation values. However, such a ground state 
is not able to explain the observed strong experimental decreaseEm of the sublattice magnetization with increasing 
hole doping. 

In this section we want to determine the ground-state energy and the staggered (sublattice) magnetization of 
the doped system. In the next section we shall show how spin-wave energies can be obtained from a ground-state 
calculation. The results demonstrate that this new static approach leads to similar results as obtained from the 
calculation of dynamical quantities (e.g., correlation functions). 

In the Hamiltonian we introduce an additional field Ba parallel to the z-axis which couples to the staggered 
magnetization. We use the following decomposition: 

H = H + E u 

Hq — H I sing ~t~ H Zeeman 

= jJ2(S?S*-^) + J(N- 2M) + 9J ^ B B A (- £ Sf + £ 5?), (26) 

<ij> i j 

Ei = H t + H ± 

= -t E ( £ t ^ + & tc ) + \ E ( s i s t + s + s i ) ■ 

<ij>,<J <ij> 

The ground state |</>o) of the unperturbed Hamiltonian Eg is an antiferromagnetically ordered Neel state with M 
holes. The holes have fixed momenta k m and are located on the sublattice a rn (a m =|, |) 



G 



\4>o) — Cfc lCT1 . . . Ck M<TM \4>Neel) 

= jwwm fi ( £ eikmR<m ( 2? ) 



The staggered magnetization at zero temperature can be obtained from the ground state energy i? by 



f) F 1 



where \tpo) denotes the exact ground state. 

In the wave operator fi we have to include both parts of Hi , the spin- flip term H± and the hopping term Ht . Spin 
fluctuations generated by H± can be described by pairs of spin waves (S~^S* ^) with 

S~ t = -i= Ye^S-, (29) 
S+ = Wst . 

defining creation operators of magnons on— both sublattices. The momenta q have to be taken from the magnetic 
Brillouin zone. Recently, we have shownE^I that the dynamics of the undoped system can well be described by 
a wave operator containing the magnon creation operators (p9[). This represents an extension of the Bartkowski 
wavefunctiont3 used in the previous section. We include spin-flip pairs not only on nearest neighbor sites but on 
sites which are arbitrarily far away from each other. To treat hole motion processes induced by H t we use the path 
operators A n from Sec. Ill without the projectors Pi 3 -. Both types of excitations are included in the wave operator fi. 



fi = exp ( " q (StfS+^y + X nA n I (30) 
\ q n / 

Note that the path operators A n commute with the spin- wave creation operators S'q (p9|), because they all contain 



only spin-flip operators destroying the Neel order. The dot ■ in the first term of (30) indicates that the quantity inside 
(...)' has to be treated as a single entity in the cumulant formation. 

The set of equations for the ground-state energy and the coefficients v q and A„ (analogous to (||)) reads 

E = (H'\Q) - ^Ng Jf i B B A (l-S), 

0= ((S-^St^Y \H' n), (31) 
0= (A n \H'Q) 

with fi given by (|3^). The cumulant expectation values have to be taken with respect to the unperturbed ground state 
@. We have shifted the energy zero level by setting H' = H + gjfj, B B A (l - S)N/2, This leads to (^ol-ff^o) = for 
M holes in the system which are not on nearest neighbor sites. After evaluating the cumulants and expanding the 
exponential we again obtain terms containing products of path operators which can again be replaced by new path 
operators. As discussed above we take care of path concatenations by replacing the set of coefficients {A„} by new 
coefficients which have to be calculated. 

To handle (|3lj) we neglect all terms which depend more strongly than linear on the hole concentration 6. In 
particular, we do not consider hole-hole interactions which are of order S 2 . This is certainly a good approximation 
for weak doping. For the energy we find 

(H\Q)= ( fl± $> q S ^r S -J + & ( + < 
q 

1^ E^ 5 qT^-q|l^)( 1 -' 5 ) 2 



q 

^Neel I 



M 



£ ( (ii (m\H t Ai\m) (1 - 6) + fa (m\H ± A 2 \m) (1 - 5) 2 ) (32) 
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where we have introduced the notation |m) as abbreviation for an one-hole state: 



\m) = c km(rfn \(t> N 6ei) ■ (33) 

The brackets (...) in the first line of j3^ ) denote expectation values with the ground state \(j>o) of Hq. In the second 
equation of (f32|) we have assumed that the holes move independently. When calculating expectation values like 
(<j>o\H±A 2 \(f>o) holes which do not couple to H± and A 2 give rise to a prefactor (1 — S) for each site where the spin 
operators from both H± and A2 act. (1 — S) describes the probability for finding a spin at such a site. For consistency 
we explicitely write down all factors (1 — S) although the actual calculation is only valid up to first order in 6. The 
energy (|32] ) consists of a spin- wave part and a hole part. The latter is proportional to the hole concentration 5 (besides 
the factors (1 — 6) mentioned above). 

For the coefficients ^ q we obtain from (|l])2 a se t of integral equations for ^ q when j3 n is fixed. Eq. ((3^)3 is a set of 
coupled non-linear equations for [3 n with fixed f q . Analogous to the preceeding section it can be transformed into a 
generalized eigenvalue problem. To be short, here we only state the integral equation for the ^ q : 

= (<f>Neel\ (^ qT ^ qi )' + H±\4>Ntel) (1 - 6) 2 

+ (<t>Neel\ (^qT^iqJ.)'" 1 " {Hl sing (l - 5) + H Z eeman) (S^S^^)' \<t>Neel) (1 - S) 2 
+ (4>mel\(S- 1 St <ii y + H ± ( S q 1T S -\u)') 2 l^eeO (1 - 



2! 

qi 



M 



+ E ( m \( s ^ s -qi)' + ( H t^ A i + {H Ismg (l - 5) + H Zeeman )p 2 A 2 ) \m) (1 - Sf 



m—1 
M 



+ E Htti)' + (¥i^ + H ± fcA 2 (l-5)) ^^ qi (5 qiT ^ qu )-|m)(l-<5) 3 

m—1 qi 

M 

+ E ( m K S ^ S -^y + (Hls lng (l -6)+ Hzeeman) ^ (Pi Al) 2 |m) (1 - 8) 2 
m—1 

M ( 1 \ 

+ E (rn\(S^S+^y+(H t f3 1 A 1 p 2 A 2 + H^-(p 2 A 2 ) 2 (l-S)j \m)(l-6f 

m=l > ' ' 

M 

+ E (™l(^ q T^qi)' + ^^(/ 3 i^i) 2 E^^q" 1 T 5 -qu)»( 1 -^ 4 • ( 34 ) 

m—1 qi 

The brackets (4>Neel\-"\<pN6el) an d (m|...|m) denote cumulant expectation values with \<pNeei) & n d \m), respectively. 
To simplify the evaluation of (34) we use linear spin- wave approximation and assume independent hole motion. 
Furthermore we cut the eigenvalue problem for the coefficients j3 n after the third variable, i.e., we include only paths 
of lengths 0, 1, and 2. Note that the essential process for the hole motion is taken into account by this approximation, 
i.e., a hole can hop twice via H t , then the spin defects are removed via H±. The mechanism of hole-spin coupling 
is also included: A hole hops via H t , the arising spin defect can be considered as starting point of a spin wave 
propagating via H±. With linear spin-wave approximation ( |34j ) reduces to a quadratic equation for each i/ q . 

To find a solution of ( |3l| ) (with fixed hole momenta) we have to proceed iteratively: starting from some fixed values 
for Pi and p 2 we first calculate the z/ q for each q. The obtained values for ^ q are then inserted into the eigenvalue 
problem for the coefficients p n . The solution leads to new values of p n for the next step. With good initial values for 
Pi and p 2 the iteration converges after a few steps. 

Taking the derivative of Eq with respect to the external field Ba according to (|2^) we can evaluate the doping 
dependend staggered (sublattice) magnetization. Results for small hole concentrations and different values oit/J are 
shown in Fig. 4. The magnetization for the undoped case becomes M e ff/M e ff,Neel ~ 0.606, i.e., it is the same as 
in linear spin-wave theoryH. This was expected since our approximations are equivalent to those of linear spin-wave 
theory. The magnetization decreases with increasing hole concentration <5. At a critical doping level 5 C1 it vanishes 
indicating a magnetic phase transition to a paramagnetic state. For instance, for t/J — 5 the cricital concentration 
S Cl is approximately 3.1 % which is in good agreement with experimental data. 

With the coefficients f q and p n we are able to compute static and dynamical properties of the system. In the next 
section we want to show how to determine spin-wave energies within our method. 
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V. SPIN-WAVE ENERGIES 



The ground state calculation in the preceeding section allows for the derivation of spin- wave energies. The basic 
idea is to calculate the energy of a state which is the ground state of the system with certain fixed boundary 
conditions. Such boundary parameters are conserved quantities like the electron number, the total magnetization 
or the total momentum. The ground state determined above has (N-M) electrons, zero total momentum and zero 
total magnetization (if we have the same number of holes on each sublattice). 

To find an expression for the spin-wave energy we now start from an unperturbed state \<fio) which contains one 
additional spin excitation 

I'M = S K ^C kllJ1 . . . Ck MCTM \<f>N6el) ■ (35) 
Since the total magnetization M tot is a conserved quantity, 

[M toU H]_ = 0, [M tot , = 0, (36) 

M tot = g., » B \J2 S i +Y1 S i ) ' ( 37 ) 



(with f2 given by the exponential form (|3C| ) ) the perturbed state has a net magnetization of gjfiB in ^-direction and 
is therefore orthogonal to the ground state calculated in the preceeding section. A straightforward calculation based 
on the wave operator ([}0]) leads to the same equations for the coefficients z^ q and [3 n . This is also due to our use of the 
spin- wave-like approximations where the spin- wave modes are treated as independent excitations. For the energy of 
the system we find an expression which differs from the ground state energy Eq of eq. ( |32| ) by (7i=l): 

AE = cjk = 2 J (1 - S) + ^y)- (1 - 5 f 

Mi 



Pi ^2{m\^H t A 1 S^\m)(l-S) 

m— 1 
M 

P* Yl MS^H ± A 2 S K] \m)(l-S) 2 . (38) 



Mj = M/2 is the number of holes on the "f-sublattice. The first m-sum runs over all holes with a m =f . Note that 



a spin excitation in |0o) (35) would lead to the same result because the hole momentum distribution in the 
ground state should be symmetric with respect to the sublattices. The quantities i/k and /3„ depending on the hole 
concentration S can be taken from the calculation in the last section. J(q) denotes the Fourier-transformed exchange 
coupling defined by 

J(q) = J 2o7 (q) = J^V qA = 2 J (cos q^ +cosq y ). (39) 

A 

For the undoped case (5 — 0) the expression (|38|) exactly reproduces the result of linear spin-wave theory: 



ujk = 2Jv/l- 7 (K) 2 . (40) 



The spin-wave energies for different hole concentrations calculated from (|38|) are shown in Figs. 5 and 6. 

For an analytic discussion of the K-dependence of the spin-wave energy we calculate z/k for fixed path coefficients 
(3\ and (3% using approximations which are again equivalent to those of linear spin-wave theory. (Note that the results 
must be inserted iteratively into the eigenvalue problem for the (3 n to get a self-consistent solution as discussed in the 
last section). Eq. ([m|) is a quadratic equation for 2/ q . Its solution can be inserted into (^). With 

M 

tA( q ):= ^M(S+Sr qi )-iWS qT ^ qi )-|m) (41) 

m— 1 

we find 







„ K - 2J (1 - V 1 - 7 (K)* - + (ft dffl '/ J)Z . (42, 

The main contribution to the renormalization of the spin-wave energy is determined by fa and is therefore due to 
the coupling between a path of length 1 and a spin wave. Note that terms with fa cancel each other. The next 
process renormalizing wk is expected to arise from fa which we neglect here. However, fa is small compared to fa, 
see Sec. III. From this fact we expect that the error introduced by including only paths of length 0,1, and 2 in the 
calculation of the spin-hole coupling is rather small. 

To perform the expectation value in A(K) we need to know the momentum distribution of the hole quasiparticlcs. 
It is well established that the minimum energy of the quasiparticle dispersion is located at (±7r/2, ±7r/2). As a simple 
approximation we neglect the effect of the Fermi surface and assume that the momenta of all holes are (±7r/2, ±7r/2), 
i.e., all holes are in the centre of the hole pockets. In the following we are only interested in small momenta K. The 
first non-trivial order of K leads to 

I(K) K 2 

7 (K) = ^-L « 1 - — , A(K) « 4<5K 2 . (43) 
J z Q 4 

Inserting this into (E^) we end up with the following approximation for the spin-wave energy at small momenta K: 



uj k = V2J |K| (1 - S) y/l - 46 fa t/J . (44) 

The K-dependence remains linear even upon doping. For K — > the spin- wave energy wk goes to zero in accordance 
with Goldstone's theorem. The spin-wave velocity v = o>k/|K| (|K| — > 0) decreases with doping. At a critical 
hole density 5 C2 the spin-wave velocity vanishes. If we assume fa to be independent of the hole concentration and 
approximate 1 — S ~ 1 we obtain 



v = v a/I - S/8 C2 , (45) 
J 



Sc2 ~~ Ufa' 

Vq denotes the spin-wave velocity in the undoped antiferromagnet. Using realistic values t/J— 5, rSi = 1.25 we find 
<5 C2 =4.0%. The doping dependence of the spin-wave velocity is in good agreement with experiments!! as well as with 
other theoretical resultsOCa. Fig. 7 shows the calculated doping dependence of v for t/J = 5. Another important 
result to be seen in Fig. 5 is that the K-range, where we find a linear K-dependence of uk, becomes smaller with 
increasing hole concentration. This is also in agreement with experimental results, see e.g. Rossat-Mignodlil. Note 
that the critical concentration S Cl where the magnetization vanishes is somewhat smaller than 8 C2 evaluated here. As 
already mentioned for t/J = 5 we have S Cl « 3.1% whereas 6 C2 « 4.0%. This is understandable since the softening of 
spin-wave excitations should first lead to the vanishing of the staggered magnetization. 

The cumulant formalism provides an alternative way to derive the expression (^) for the spin-wave energy. Intro- 
ducing a transverse wave vector dependent field in the Hamiltonian one obtains the transverse static susceptibility by 
taking the second derivative of the ground state energy Eq with respect to the field. The results for the static suscepti- 
bility allow for the calculation of spin-wave energies. The necessary link is given by the dynamical spin susceptibility. 
The calculation leads exactly tevthe same expression ([38]) for the spin- wave energy (!). The detailed derivation has 
been published in a recent papeiH. The static susceptibility for a staggered field shows a K _2 -divergency for all hole 
concentrations. This means that the calculation presented in Sec. IV always describes a state with antiferromagnetic 
long-range order. 



VI. CONCLUSION 



The aim of this work was to study the hole and spin dynamics of weakly doped antiferromagnets described by 
the t-J model. Instead of using dynamical correlation functions our approach is based on the calculation of the 
ground-state energy and contains a static view of the system. 

In the first part we have rederived the quasiparticle dispersion for one hole generated in the ground state at half- 
filling (Fig. 1,2). The energy minima are found at (±7r/2, ±7t/2) in agreement with other analytical and numerical 
results. By neglecting hole-hole interaction the calculation can be extended to finite doping by use of a rigid-band 
approximation. For the Fermi surface (Fig. 3) we obtain the expected picture of hole-pockets around (±7r/2, ±7r/2). 

The weakly doped system has been subject of the second part of this paper. We have presented results for the 
staggered magnetization and for spin-wave energies up to first order in 6. The staggered magnetization (Fig. 4) 
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decreases with doping due to spin-hole interactions. The magnetization becomes zero at a hole concentration S Cl 
which indicates the disappearance of antiferromagnetic long-range order. The spin- wave energy is found to be also 
strongly renormalized with hole doping (Fig. 5 and 6). Both effects are coupled: Due to decreasing spin-wave energies 
more long-wavelength spin fluctuations are mixed into the ground state which causes the loss of magnetization. The 
spin-wave velocity has a square-root concentration dependence (Fig. 7) and vanishes for a critical hole concentration 
S C2 given by ((f5[). Note that we find S Cl < S C2 (for instance, S Cl = 3.1%, S C2 = 4.0% for t/J = 5), i,e., the magnetization 
vanishes before the spin-wave velocity becomes zero. This feature was observed in experiments!!!. 

The strong renormalization of the spin-wave energies due the presence of holes arises from the interaction of spin 
waves with spin fluctuations created by hole hopping. This can be interpreted as the decay of spin waves into particle- 
hole pairs. To first order in 6 only the first path coefficient /3i contributes to the renormalized spin- wave energy ([42]). 
The destruction of the antiferromagnetic state therefore results from the creation and annihilation of additional spin 
waves by moving holes. The strength of this coupling increases with t/J because the energy gain by hole hopping is 
of order t while the energy loss due to a flipped spin is of order J. Realistic values for t/J are between 3 and 5 which 
explains the high effectiveness of a small hole concentration for the destruction of antiferromagnetism. 

In principle it should also be possible to discuss the influence of spin dynamics on the quasihole properties. However, 
the approximation of truncating the eigenvalue problem for the path coefficients after the third variable (A?) causes 
the loss of some features of the one-hole dispersion relation. This approximation leads to the fact that all momenta 
on the boundary of the magnetic Brillouin zone correspond to the same energy, i.e., in this case the minima are no 
longer obtained at (±7r/2, ±7r/2). It would be necessary at least to include also paths with a length of 4. This would, 
however, lead to a more extensive calculation. 

The present calculations are only valid within the antiferromagnetic phase of the t-J model at small doping con- 
centrations. Our starting point is the Neel state which has long-range order. Spin fluctuations are included by a 
perturbational method based on cumulants. Note that the linear spin-wave approximation employed for calculating 
the expectation values becomes questionable for vanishing sublattice magnetization. Therefore, the actual calcula- 
tion can not provide a reliable description of the system in the vicinity of the magnetic phase transition (S « S Cl ). 
Especially properties of the critical point are not accessible. 

In the present approach, even for higher doping the ground state has always a broken rotational symmetry. The 
divergency of the staggered susceptibility for small momentat3 ~ 1 /K 2 shows that we always describe a system with 
antiferromagnetic long-range order. To extend our approach to systems with only short-range magnetic order will 
be the subject of future research. 



APPENDIX A: EVALUATION OF CUMULANTS 



In this appendix we show how to evaluate cumulants with an exponential ansatz for the wave operator f2, compare 
(^). The following identities are useful for transforming the set of non-linear equations ((2^)2 into an eigenvalue 
problem. 

(Ae s ) c = (Ae s ) , 
(ABc s ) c = (ABe s ) - (Ae s ) (Be s ) , 
(Ae s B) c = (Ac s B) - (Ae s ) (e s B) , 
(ABCc s ) c = (ABCe s ) - (ABe s ) (Ce s ) - (ACe s ) (Be s ) 
-(BCe s ) (Ae s ) + 2(Ae s ) (Be s ) (Ce s ) , 

(Al) 

These relations hold for operators with (S ) — OVfc > 0, (A) = 0, (B) = etc. They can be proven straighforwardly 
by expanding the exponentials, explicitely evaluating the cumulants and then recollecting all terms. To show this for 
the second identity we start from the definition of cumulant expectation valuesEj for a product of arbitrary operators 
At 

(0in^> c = (nf^r) )M<f>\u eXiAi w\^i- ( a2 ) 



We find 



(ABS k ) c = ±JL(±) ln(e*W s ) | A=K=J)=0 



dX dn \dr] 
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8 \ k (ABe" s )(e" s ) - (Ae r > s )(Be , ' s ) 



ch] 



1 17=0 



= {ABS k ) - (*) (AS^iBS' 



Here we have used (S , ' c ) = 0. Inserting this transformation in the original expression leads to 



(ABe s ) c = (AB) C + J ±.(ABS k ) c 
fe=i 

oo oo 1 

= (ABe s ) — ^2 ^2 — (AS i )(BS j ) 



i=0 j=0 



(A3) 



(A4) 



In the third equation we have replaced the sum over k by a sum over j with j = k — i. Eq. (A4) is the desired result. 
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FIG. 1. Energy dispersion of the lowest hole excitation for different values of t/J. The zero energy was set to the center of 
mass of the band. 
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FIG. 2. Total bandwidth of the hole dispersion shown in Fig.l as function of t/J 
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FIG. 3. Hole Fermi surface obtained from the dispersion relation (Fig. 1) in a rigid band approximation for hole concentrations 
5=5% and 15%, t/J=5. 
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